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This  dissertation  presents  the  results  of  a study  on  the  nonlinear  stochastic  analysis  of 
random  waves  in  finite  water  depth,  in  particular,  (1)  clarification  of  nonlinearity  in  wave 
spectra  and  (2)  development  of  a probability  density  function  in  closed  form  applicable  to 

wave  profiles. 

In  order  to  clarify  the  nonlinear  characteristics  associated  with  wave-wave  interactions, 
which  are  particularly  pronounced  for  waves  in  finite  water  depth,  a method  is  developed  to 
examine  those  frequencies  having  nonlinear  energy  components  and  their  magnitudes.  The 
nonlinear  components  of  the  spectral  density  at  a specified  frequency  are  considered  to  be 
the  accumulation  of  nonlinear  interactions  associated  with  various  pairs  of  two  frequency 
components.  Here,  interaction  includes  not  only  the  sum  but  also  the  difference  of  two 
frequencies.  The  separation  of  the  nonlinear  component  of  the  energy  density  of  a spectrum 

is  achieved  by  applying  the  concept  of  the  bicoherence  spectrum. 

The  separation  procedure  is  applied  to  wave  spectra  computed  from  data  obtained  in 
the  ARSLOE  Project  during  storm.  The  results  demonstrate  that  nonlinear  components  are 
present  at  low  and  high  frequencies  but  no  nonlinear  components  exist  in  the  neighborhood 
of  the  frequency  where  the  spectrum  peaks  and  that  the  ratio  of  nonlinear  energy  to  the 
total  energy  increases  significantly  with  a decrease  in  water  depth. 


IX 


In  order  to  evaluate  the  statistical  properties  of  nonlinear  waves,  a probability  density 
function  applicable  to  non-Gaussian  random  waves  is  developed  in  closed  form.  Currently, 
no  probability  density  function  in  closed  form  representing  non-Gaussian  random  processes 
is  available. 

In  the  derivation  of  the  probability  density  function,  the  concept  of  Kac-Siegert’s 
method  developed  for  nonlinear  mechanical  system  is  applied.  That  is,  the  Kac-Siegert 
formula  is  asymptotically  expressed  in  terms  of  a random  variable  that  obeys  a normal 
distribution  with  parameters  evaluated  from  information  on  cumulants  of  the  wave  record. 
Then,  by  applying  the  transformation  technique  of  random  variables,  the  desired  probabil- 
ity density  function  is  developed  in  closed  form.  Comparisons  between  the  newly  developed 
probability  density  function  and  the  histograms  constructed  from  wave  records  at  various 
water  depths  obtained  in  the  ARSLOE  Project  during  storm  show  excellent  agreement. 


X 


CHAPTER  1 
INTRODUCTION 


The  characteristics  of  waves  in  finite  water  depth  are  significantly  different  from  those 
of  waves  in  deep  water.  Visibly,  waves  in  finite  water  depth  show  a trend  of  skewness  in 
their  profile,  and  the  degree  of  skewness  is  pronounced  in  severe  seas  in  shallow  water. 
This  is  attributed  to  a nonlinear  wave- wave  interaction  (namely,  energy  transfer)  between 
component  waves. 

Although  many  studies  on  the  nonlinear  characteristics  of  ocean  waves  have  been  carried 
out  to  date,  little  research  has  been  made  on  nonlinear  stochastic  properties  of  waves.  Those 
studies  are  primarily  based  on  the  perturbation  expansion  for  a random  sea  of  constant  finite 
depth  such  as  the  weakly  nonlinear  random  waves  derived  from  the  Stokes’  expansion  of 
governing  equation  using  potential  theory  (Tick  1959,  1961,  Phillips  1960,  1961,  Hasselmann 
1962.  1963a,  1963b,  Longuet- Higgins  1963).  The  weakly  nonlinear  wave  theory  has  provided 
significant  insight  on  nonlinear  wave- wave  interactions;  however,  the  theory  cannot  account 
for  steep  waves  or  waves  in  shallow  water  that  show  strong  nonlinearity. 

It  is  of  considerable  interest  to  understand  nonlinearity  involved  in  the  spectra  of  waves 
in  finite  water  depth.  Several  studies  on  the  quantitative  estimation  of  nonlinear  spectral 
energy'  components  of  ocean  waves  have  been  conducted  (Masuda  et  al.  1979,  Hudspeth  and 
Chen  1979,  etc.).  However,  almost  all  these  studies  have  been  based  on  weakly  nonlinear 
random  waves,  which  are  obviously  insufficient  to  explain  nonlinearity  in  the  spectra  of 

waves  in  finite  water  depth. 

Another  important  and  interesting  study  on  nonUnear  characteristics  of  waves  in  finite 
water  depth  is  the  probability  distribution  of  the  wave  profile.  Currently  available  proba- 
bility density  functions  representing  non- Gaussian  characteristics  of  wave  profiles  m finite 
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water  depth  have  drawbacks.  For  example,  the  probability  density  functions  given  in  series 
form  (Gram-Charilier  prob.  distribution)  have  negative  values  at  some  part  of  the  distri- 
bution. The  probability  density  functions  developed  based  on  the  weakly  nonlinear  wave 
theory  cannot  be  applied  to  waves  having  strong  non-Gaussian  characteristics  (Tayfun  1980, 
Huang  et  al.  1983,  Langley  1987). 

This  study  adresses  nonlinear  characteristics  of  waves  in  finite  water  depth,  in  particu- 
lar, (1)  clarification  of  the  nonlinearity  in  wave  spectra  and  (2)  development  of  a probability 
density  function  in  closed  form  that  is  applicable  to  surface  elevation  of  waves  in  finite  water 

depth. 

In  order  to  clarify  the  nonlinearity  in  wave  spectra,  a bispectral  analysis  technique  is 
applied  for  quantitative  evaluation  of  nonlinear  wave-wave  interactions.  When  various  spec- 
tral components  interact  with  one  another  due  to  nonlinearities,  the  resulting  components 
are  phase  coherent  with  the  primary  components.  Therefore,  those  components  are  not 
considered  to  be  statistically  independent.  The  bispectrum,  which  measures  the  statistical 
dependence  of  three  spectral  components  whose  sum  frequncy  is  zero,  may  therefore  be  a 
useful  tool  to  investigate  the  nonlinearities  of  random  waves  in  finite  water  depth.  By  apply- 
ing the  concept  of  a bicoherence  spectrum,  a method  to  separate  the  nonlinear  component 

of  the  energy  density  of  spectrum  will  be  developed. 

In  the  derivation  of  a probability  density  function  applicable  to  nonlinear  wave  profiles, 
the  concept  of  the  Kac-Siegert  method  (Appendix  A)  for  evaluating  the  response  of  a 
nonhnear  mechanical  system  is  applied.  The  background  of  the  application  is  as  follows: 

It  is  commonly  known  that  the  response  of  a nonlinear  system  is  generally  a non- 
Gaussian  random  process,  even  if  the  input  of  the  system  is  a Gaussian  random  process.  The 
response  of  a nonlinear  system  is  often  presented  by  Volterra’s  stochastic  series  expansion, 
and  Kac-Siegert  derived  the  probability  distribution  of  the  response  with  knowledge  of  the 
second-order  frequency  response  function. 

Now,  let  us  consider  wave  profiles  in  the  ocean.  The  results  of  statistical  analysis  of 
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wa.v6s  obtctincd  during  the  Atlantic  occcin  R,6motG  Sensing  Lund-Ocenn  Experiment  (AR 
SLOE)  Project  (Ochi  et  al.  1983)  have  shown  that  the  profile  of  offshore  waves  follows  the 
normal  probability  distribution,  while  the  wave  profile  near  shoreline  obeys  the  Gram- 
Charlier  series  (non-Gaussian)  distribution.  This  implies  that  wave  characteristics  are 
changing  from  a Gaussian  random  process  to  a non-Gaussian  random  process  with  decreas- 
ing water  depth.  This  situation  resembles  the  response  of  a nonlinear  mechanical  system 
which  is  a non-Gaussian  random  process  whose  skewness  increases  with  increase  in  degree 
of  nonlinearity.  The  concept  of  Kac-Siegert’s  method  which  yields  the  probability  density 
function  of  a nonlinear  response  may  therefore  be  applied  to  the  problem  at  issue.  However, 
the  second-order  frequency  response  function  necessary  to  apply  the  method  is  unknown 
for  nonlinear  waves.  Hence,  we  will  develop  an  analytical  method  to  derive  Kac-Siegert  s 
solution  from  only  a knowledge  of  nonlinear  wave  data  through  separating  wave  spectra 
into  linear  and  nonlinear  components. 

The  paper  consists  of  three  chapters  followed  by  conclusions  and  appendices.  Chapter 
2 reviews  currently  available  literature  on  nonlinear  components  in  wave  spectra  as  well  as 
the  probability  distributions  applicable  to  surface  elevation  of  nonlinear  random  waves. 

Chapter  3 presents  a method  to  identify  second-order  spectrum  representing  the  non- 
linear properties  in  energy  spectrum  computed  from  mecisured  wave  data.  The  method  is 
applied  to  field  data  obtained  in  the  ARSLOE  Project  during  storm. 

Chapter  4 presents  the  derivation  of  the  probability  density  function  in  closed  form 
applicable  to  random  waves  in  finite  water  depth,  which  is  considered  to  be  a non-Gaussian 
random  process.  By  applying  the  concept  of  Kac-Siegert’s  method,  a new  probability  den- 
sity function  is  derived  in  closed  form  and  is  compared  with  histograms  constructed  from 
measured  wave  data. 

Chapter  5 presents  the  findings  and  conclusions  obtained  from  the  present  study. 

The  appendices  outUne  (A)  Kac-Siegert’s  method,  (B)  ARSLOE  wave  data,  (C)  the 
computational  procedure  of  bispectrum,  and  (D)  inverse  presentation  oi  Y = U -[■  . 


CHAPTER  2 
LITERATURE  REVIEW 


This  chapter  reviews  literature  on  nonlinear  components  in  wave  spectra  and  the  prob- 
ability distributions  applicable  to  nonlinear  random  waves,  which  usually  can  be  considered 
as  a non- Gaussian  random  process. 

2.1  Nonlinear  Components  in  Wave  Spectra 

The  perturbation  expansion  (or  Stokes’  expansion)  method  to  the  governing  equations 
of  inviscid  fluid  and  incompressible,  irrotational  flow  has  been  extended  to  random  sea 
waves  by  Tick  (1959,  1961),  Phillips  (1960,  1961),  Hasselmann  (1962,  1963a),  and  Longuet- 
Higgins  (1963).  Hereafter,  the  solution  obtained  from  the  Stokes’  expansion  method  for 
random  waves  in  finite  water  depth  will  be  denoted  as  the  weakly  nonlinear  solution. 

Tick’s  derivation  of  the  weakly  nonlinear  solution  extended  to  the  second  order,  and 
he  presented  the  second-order  interaction  kernel  for  random  waves  in  finite  water  depth 
(1959,  1961).  Hasselmann  (1962)  carried  out  the  perturbation  to  the  third  order  in  spectral 
energies,  and  found  that  the  spectral  energy,  could  be  represented  in  a perturbation 
series  as  follows: 

E = J?2  + ^^4  + E^  + * * • (2*1) 

where 

E2  = linear  spectrum,  spectral  energy  due  to  linear  components 
E4  = second-order  spectrum 
Ee  = third-order  spectrum 

The  perturbation  terms  of  odd  order  all  vanish  for  a Gaussian  zero  mean  varible. 

E3  = Es=--  = 0 (2.2) 
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Since  the  spectrum  obtained  from  measured  wave  data  contains  both  linear  and  nonlin- 
ear higher  order  spectral  energy  components,  it  is  often  necessary  to  separate  the  spectral 
energy  of  linear  components  from  the  measured  spectrum.  Hereafter,  the  spectrum  obtained 

from  measured  data  will  be  denoted  as  the  measured  spectrum. 

The  second-order  spectrum  52(w)(=  E4)  is  given  by  the  following  convolution  integral 

(Tick  1961,  for  example). 


/OO 

jiir(a>fc,u;)p5i(a;  - ujk)Si{uk)duJk 

-00 

where 


SiM 

U>jt,  w 


E2  = linear  spectrum 
second-order  interaction  kernel 
radian  frequency 


Many  studies  have  been  conducted  to  separate  the  second-order  spectrum  from  the 
measured  spectrum.  As  shown  in  Eq.  2.3,  the  kernel  function  and  linear  spectrum  should  be 
given  to  compute  the  second-order  spectrum.  Most  of  the  studies  reviewed  here  employed 
the  second-order  interaction  kernel  obtained  from  the  weakly  nonlinear  solution.  Since 
the  linear  spectrum  5i(u;)  in  Eq.  2.3  is  unknown  in  advance,  each  method  has  a different 

approach  in  assuming  the  first-order  spectrum. 

The  first  attempt  to  compute  the  second-order  spectrum  was  made  by  Tick  (1959, 1961). 
In  the  derivation  of  the  second-order  weakly  nonlinear  solution,  he  computed  the  second- 
order  spectrum  by  assuming  that  the  linear  spectrum  was  given.  Hamada  (1965)  made 
a correction  to  Tick’s  solution.  To  obtain  the  second-order  spectrum  from  the  measured 
spectrum  in  the  wave  tank,  he  assumed  the  approximate  shape  of  the  linear  spectrum 
by  excluding  the  part  that  appeared  to  be  the  obvious  nonlinear  part  from  the  measured 
spectrum.  In  order  to  compute  the  second-order  spectrum,  Hudspeth  and  Chen  (1979) 
employed  both  theoretical  spectra  and  a smoothed  measured  spectrum  estimated  from  wave 
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data  obtained  during  Hurricane  Carla.  Masuda  et  al.  (1979)  developed  the  iterative  method 
for  the  separation  of  the  second-order  spectrum  from  the  measured  spectrum.  They  assumed 
that  the  measured  spectrum  consisted  of  linear  and  second-order  spectrum.  Using  the 
measured  spectrum  as  an  initial  estimation  for  the  linear  spectrum,  they  could  separate  the 

second-order  spectrum  within  a few  iterations. 

The  methods  for  obtaining  the  second-order  spectrum  reviewed  above  are  all  based  on 
the  second-order  interaction  kernel  of  weakly  nonlinear  solution.  Therefore,  as  stated  in 
Chapter  1,  these  methods  are  not  applicable  to  random  waves  having  strong  non-Gaussian 

characteristics. 

On  the  other  hand,  Kim  and  Power  (1979)  utiUzed  bispectral  analysis  to  investigate 
the  nonlinear  wave-wave  interactions  of  coherent  waves  in  plasma  fluctuation  data.  In 
particular,  Kim  and  Power  successfully  showed  that  the  bicoherence  spectrum  might  be  used 
to  discriminate  between  nonlinearly  coupled  waves  and  linear  free  waves.  Furthermore,  the 
bicoherence  squared  was  shown  to  measure  the  fraction  of  wave  energy  due  to  the  second- 
order  wave  coupling  in  a measured  spectrum.  However,  they  considered  only  the  case  of 
single  three  wave  coupling  and  did  not  apply  the  technique  to  the  separation  of  the  second- 

order  spectrum  from  the  measured  spectrum. 

The  bispectral  analysis  technique  by  Kim  and  Power  (1979)  will  be  reviewed  in  detail 
in  Chapter  3 because  the  technique  can  be  used  for  identifying  the  nonlinear  components 

of  random  waves  having  strong  non-Gaussian  characteristics. 

2.2  Probability  Distributions  Applicable  to  Nonlinear  Random  Waves 

One  of  the  early  attempts  to  represent  the  probabilistic  distributions  of  the  surface 
elevation  of  nonlinear  random  waves  was  made  by  Kinsman  (1960)  by  means  of  a Gram- 
Charlier  probability  density  function.  The  Gram-Charlier  probability  density  function  was 
developed  by  applying  the  concept  of  orthogonal  polynomials  with  respect  to  the  probability 
density  function.  Since  Hermite  polynomials  are  orthogonal  with  respect  to  the  Gaussian 
(normal)  distribution,  the  probability  density  function  can  be  expressed  as  a function  of 
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Gaussian  distribution  and  a series  of  Hermite  polynomials.  The  Gram-Charlier  series  of 
Type  A probability  density  function  for  a standardized  random  variable  is  a^  follows: 


f{z)  = 


+ 


1 


, , rr  / \ 3 IT  / \ - lOma 

1 + —Hz{z)  + — 7i — + 


3! 


4! 


5! 


H^{z)  (2.4) 


me  — 15m4  + 30 

6! 


Heiz)  + 


where 


z — standardized  random  variable 
rrij  = j-th  moment  of  the  standardized  random  variable 
Hn{z)  = Hermite  polynomials  of  degree  n 

Significant  advances  in  theoretical  analysis  of  nonlinear  random  waves  were  made  in 
the  late  1950s  and  early  1960s  by  Tick  (1959,  1961),  Phillips  (1960,  1961),  Hasselmann 
(1962,  1963a,  1963b),  and  Longuet-Higgins  (1963).  Longuet-Hig^ns  (1963)  specifically 
investigated  the  effects  of  nonlinearities  on  the  statistical  distribution  of  random  waves.  He 
derived  a series  probability  density  function  through  the  cumulant-generating  function.  The 
distribution  was  expressed  in  terms  of  cumulants  as  foUows: 


where 

cr^  = variance  of  y 

3 

A3  = ks/k^  = skewness 
A4  = k^lk\  = kurtosis  — 3 

As  = ksikl 
k * — 


cumulants 
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This  expression  proved  to  be  in  the  same  form  as  the  Edgeworth  series  distribution 
(1905)  and  also  the  modified  form  of  the  Gram-Charlier  probability  density  function.  There- 
fore, the  Gram-Charier,  Edgeworth,  and  Longuet-Higgins  probability  density  functions  were 
essentially  the  same  form,  although  their  approaches  were  entirely  different.  Detailed  com- 
parisons and  derivations  of  these  three  probability  density  functions  were  presented  by  Ochi 

(1985). 

As  is  shown  in  Eq.  2.5,  the  first  term  is  the  normal  distribution  and  the  remaining  terms 
represent  the  effects  of  non-Gaussian  characteristics.  This  probabilty  density  function, 
however,  has  a drawback.  Since  the  distribution  is  given  in  series  form,  it  usually  has 
a negative  value  at  some  part  of  the  distribution  caused  by  a finite  number  of  terms  in 
computing  the  density  function.  Ochi  and  IVang  (1984a)  applied  the  probability  density 
function  to  wave  records  measured  at  various  water  depths  during  the  growing  stage  of  a 
storm  and  concluded  that 

• the  series  is  not  uniformly  convergent;  in  other  words,  the  inclusion  of  higher  order 
terms  does  not  necessarily  give  more  accurate  results,  and 

• the  probability  distribution  that  includes  the  parameter  A3  (skewness)  only  agrees 
relatively  well  with  the  histogram,  in  general. 

Another  approach  for  deriving  the  probability  distributions  applicable  to  nonlinear 
random  waves  is  from  the  weakly  nonlinear  solution.  A second-order  random  wave  may  be 
expressed  in  the  following  form: 

y{t)  = iZef;cfce'^‘  + iZe  f]  c^c,  -H  (2.6) 

fc=l  A:=l  /=1 

where 

(j)f^  = uJkt  + €k  = phase  function 

t = time 
Uk  = 


27rA:/T  = radian  frequency 
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T = record  length  of  time  series 

£jfe  = uniformly  distributed  random  phase 

qki  is  a coefficient  associated  with  Wfc  + wj,  which  is  the  sum  of  the  interactions  between  two 
frequency  components  Uk  and  w/,  while  rjt/  is  associated  with  the  difference  of  interactions 
between  Wfc  and  wj.  These  interaction  coefficients  are  functions  of  wjt,  w/,  and  the  water 
depth.  Eq.  2.6  has  been  used  extensively  for  the  digital  simulation  of  nonlinear  random 
waves  (Hudspeth  and  Chen,  1979;  Sharma  and  Dean,  1979;  Anastasiou  et  al.,  1982;  Tuah 
and  Hudspeth,  1985)  and  for  statistical  analysis  (Tayfun,  1980;  Langley,  1987). 

Tayfun  (1980)  developed  the  probability  density  function  applicable  to  noidinear  ran- 
dom waves  observed  in  deep  water.  Assuming  the  spectrum  was  narrow-banded,  he  devel- 
oped from  Eq.  2.6  the  following  amplitude  modulated  Stokes  waves: 

y{x,  t)  = a(x,  t)  cos(xo  + f)  + ^koa^(x,  t)  cos  2(xo  + £)  + •••  (2.7) 

where 

X = horizontal  axis 
ko  = ^olg  = wave  number 
cjo  = mean  frequency 
g = gravitational  acceleration 

Tayfun  assumed  that  the  slowly  varying  amplitude  a(x,/)  and  the  phase  e follow  the 
Rayleigh  and  uniform  probability  distributions,  respectively,  and  derived  the  marginal  cu- 
mulative distribution  of  surface  elevation,  F{y),  by  applying  the  technique  of  transformation 
of  random  variables.  However,  the  cumulative  distribution  F{y)  and  the  probability  density 
function  f{y)  require  numerical  integration.  Therefore,  the  probability  density  function  is 
not  in  closed  form. 

Huang  et  al.  (1983)  also  studied  nonlinear  characteristics  of  waves  and  extended  Tay- 
fun’s  analysis  to  the  third  order  by  using  regular  Stokes  waves  in  deep  water  as  follows: 
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1 ci^k  ScL^k^ 

y{x,t)  = -a^k  + a cos  X + — cos  2%  H ^ — cos  3%  + 


(2.8) 


By  using  the  naxrow-bajid  assumption  and  by  analogy  with  Xayfun  s derivation,  they  as- 
sumed amplitude  a and  phase  x to  be  slowly  varying  functions.  Then,  by  transformation 
of  random  variables  with  successive  approximations,  the  probability  density  function  of  the 
surface  elevation  was  derived.  In  precisely  the  same  manner  as  in  the  deep  water  case, 
Huang  et  al.  (1983)  derived  the  probability  density  function  of  the  surface  elevation  in  finite 
water  depth  by  extending  the  method  to  the  Stokes  waves  of  the  third  order  in  finite  water 
depth,  which  was  given  by  Bowden  (1948)  as  follows: 


y = 


+ 


a?k 

a^k 


coth  kd  + a cos  x 


(2.9) 


+ 


2 

3a^k^ 

8 


coth  fcd  ( 1 -t- 


2sinh^  kd 


cos2x 


1 + - 


sinh^  kd  sinh”*  kd  8 sinh®  kd 


+ 


cos  3x 


where  d is  the  water  depth.  Huang  et  al.  found  that  the  probability  density  function  derived 
for  deep  water  waves  depended  on  two  parameters:  the  root-mean-square  surfax:e  elevation 
and  the  significant  slope.  For  waves  in  finite  water  depth,  an  additional  parameter,  the 


nondimensional  depth,  was  also  required. 

Langley  (1987)  presented  the  probability  distribution  of  second-order  random  waves  in 
finite-water  depth  without  the  narrow-band  assumption.  In  the  development  of  the  density 
function  of  wave  profiles,  he  first  represented  the  second-order  random  waves  of  Eq.  2.6 
in  the  following  form  by  using  the  technique  developed  by  Kac  and  Siegert  (1947)  for  the 
statistical  analysis  of  nonlinear  systems. 


yit)  = J2{f3iZj  + XjZj)  (2.10) 

i 

where  Zj  is  a standardized  normal  random  variable.  The  parameters  (3j  and  Xj  are  evaluated 
from  matrices  obtained  from  the  weakly  nonlinear  solution. 
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Since  y{t)  is  represented  in  terms  of  the  standardized  normal  random  variable  Zj,  the 
probability  density  function  of  y{t)  can  be  derived  through  the  characteristic  function  of 
Eq.  2.10,  but  the  probability  density  function  is  not  expressed  in  closed  form. 

It  should  be  noted  that  the  probability  density  functions  derived  by  Tayfun  (1980), 
Huang  et  al.  (1983),  and  Langley  (1987)  from  the  weakly  nonlinear  solution  are  generally 
not  applicable  to  steep  waves  or  shallow-water  waves.  This  is  because  the  probability 
density  function  was  derived  from  the  solution  of  the  Stokes  expansion  method,  which  is 
valid  only  for  weakly  nonlinear  and  nonresonant  waves.  For  example,  the  second-order 
solution  represented  in  Eq.  2.6  becomes  divergent  without  bound  in  shallow  water  where 
the  resonant  condition  is  easily  satisfied  for  second-order  wave-wave  interactions. 

Recently,  Cieslikiewicz  (1990)  applied  the  maximum  entropy  method  to  nonlinear  ran- 
dom waves.  The  maximum  entropy  method  has  often  been  considered  in  determining  the 
probability  density  function  of  a random  variable  from  data  in  which  information  on  the  sta- 
tistical distribution  has  been  inadequate.  The  principle  of  the  method  is  that  the  probability 
distribution  that  maximizes  the  entropy  function  subject  to  moments  constraints  has  the 
least  biased  statistical  information  (Jaynes,  1957;  Dowson  and  Wragg,  1973;  Sobczyk  and 
Trebicki,  1990).  For  given  moments  up  to  fourth  order  evaluated  from  data,  Cieslikiewicz 
derived  the  probability  density  function  numerically  by  solving  the  system  of  nonlinear  inte- 
gral equations.  Although  the  results  showed  good  agreement  with  field  data,  the  approach 
was  essentially  a data  fitting  technique.  One  should  also  note  that  the  probability  den- 
sity function  does  not  exist  when  moments  of  data  only  up  to  third  order  (variance  and 
skewness)  are  available. 

As  this  brief  review  of  the  literature  indicates,  most  studies  on  nonlinear  components 
in  wave  spectra  have  employed  the  weakly  nonlinear  solution  as  a tool  to  investigate  non- 
linearities  in  random  waves  in  finite  water  depth.  This  approach  is  valid  only  for  sbghtly 
non-Gaussian  processes.  The  bispectral  analysis  technique  has  proven  to  be  a useful  tool 
for  identifying  nonlinear  components  in  random  waves  and  will  be  used  in  this  paper  for 
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the  statistical  analysis  of  random  waves  in  finite  water  depth. 

Several  probability  density  functions  have  been  proposed  to  represent  the  non-Gaussian 
characteristic  of  a random  sea  surface.  The  Gram-Charlier,  Edgeworth,  and  Longuet- 
Higgins  probability  density  functions  are  essentially  the  same  form  and  have  negative  value 
at  some  part  of  the  distribution.  The  probability  density  functions  developed  based  on  the 
analytic  solution  of  Stokes  perturbation  method  are  generally  not  applicable  to  steep  waves 
or  shallow-water  waves.  The  probability  density  function  obtained  by  applying  the  maxi- 
mum entropy  method  is  essentially  a data  fitting  technique.  In  short,  probability  density 
functions  applicable  to  nonlinear  random  waves  having  strong  non-Gaussian  characteristics 
have  not  been  analytically  developed  to  date. 


CHAPTER  3 

EVALUATION  OF  NONLINEARITY  IN  WAVE  SPECTRA 


This  chapter  consists  of  two  sections.  The  first  section  presents  the  method  to  separate 
the  second-order  spectrum  representing  nonlinear  properties  from  measured  spectrum.  A 
method  for  evaluating  wave-wave  interaction  between  two  specified  frequencies  developed 
by  Kim  and  Power  (1979)  is  applied  to  nonlinear  interactions  associated  with  various  pairs 

of  frequency  components  in  the  wave  spectrum. 

The  second  section  presents  sample  application.  The  method  developed  is  applied  to 
field  data  obtained  during  the  growing  stage  of  a storm  in  the  ARSLOE  project.  The  results 
of  the  analysis  and  a discussion  on  nonlinear  characteristics  of  random  waves  in  different 
water  depths  are  presented. 

3.1  Separation  of  Second-order  Spectrum  from  Measured  Spectrum 

In  this  study,  the  time  history  of  a wave  profile  divides  into  linear  and  second-order 
part  as  follows: 

y(0  = yi(0  + J/2(0  (3-1) 

Let  us  first  consider  the  interaction  associated  with  two  frequency  components  wjt  and 
ui,  which  are  arbitrarily  chosen  constant  frequencies.  In  this  case,  Kim  and  Power  (1979) 
derived  the  interaction  in  the  spectral  density  at  a frequency  Um  (where  Um  is  the  sum  of 
(jjk  and  a;/)  as  follows: 


5(w,n)  = E [lYiMf 


+ 


1-B(u;jk,u;,)l' 


(3.2) 


where 


y^(cjrn)  = Fourier  transform  of  the  linear  component  yi{t)  at  the  frequency  Urn 
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Y{uJk)^  Y[ui)  = Fourier  transform  of  the  wave  profile  y{t)  at  the  frequencies  Uk  and  cj/, 

respectively 

= bispectral  density  function 

The  first  term  of  Eq.  3.2  represents  the  spectral  density  of  the  linear  component  yi{t) 
at  the  frequency  Umj  while  the  second  term  represents  that  of  the  nonlinear  component.  By 
multiplying  = E[\Y {(jJm)\^]  to  the  numerator  and  denominator,  respectively,  of  the 

second  term  of  Eq.  3.2,  we  may  write 


Siu^m)  = E[\YlM\^]  + 


E[\Y{uk)Y{u,)\^]  E[\YM\-^] 


S{0Jm) 


(3.3) 


= E [lFjr(u;„)p]  + b\uk,ui)  ■ Sium) 


where  ^ 

,(  x_f Y ^3  4', 

b{uJk,^,)  £[|r(u,„)|2]  / 

b((jJk,0Ji)  in  the  above  equation  is  called  the  bicoherence  spectrum,  and  its  squared  quantity 
indicates  the  fraction  of  the  spectral  density  considered  as  a nonlinear  component. 

In  applying  Kim  and  Power’s  method  to  the  present  problem,  we  consider  the  frequency 
u>k  and  ui  not  to  be  constant;  instead,  they  are  variables.  This  implies  that  the  nonlinear 
component  of  the  spectral  density  at  a frequency  u>m  is  equal  to  the  accumulation  of  non- 
linear interactions  a.ssociated  with  various  pairs  of  frequency  components  u>k  and  u>i  under 
the  condition  that  Uk  Y Furthermore,  we  consider  the  interactions  that  occurs 

not  only  at  the  frequency  (uk  + w/)  but  also  at  the  frequency  {uk  - <^i),  where  u>k  > u>i. 
The  latter  is  equivalent  to  the  interaction  between  u>i  and  (wjt  — a;;).  In  evaluating  the 
interaction  at  the  frequency  {uk  it  is  assumed  that  the  spectral  energy  density  at  fre- 

quencies smaller  than  the  minimum  frequency  of  the  main  energy  density  in  the  spectrum 
(u>s  itt  Figure  3.1)  is  due  to  the  nonlinear  interaction  associated  with  the  difference  between 
various  combinations  of  the  two  frequency  components  at  u>k  and  a>(.  The  spectral  density 


SPECTRAL  DENSITY  IN  - SEC 
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Figure  3.1:  Definition  of  frequency  u, 

for  frequencies  greater  than  w,  is  due  to  the  nonlinear  interaction  associated  with  the  sum 

of  various  combinations  of  two  frequency  components. 

Let  us  evaluate  the  interaction  due  to  the  sum  of  two  frequency  components.  We  may 

write  the  Fourier  transform  of  y(t)  at  the  frequency  Um  as  follows: 

v(ujm)  = Ai(uk,uji)yL(^k)yLM  (^-s) 

“♦"^1  — 

Note  that  the  summation  is  for  various  combinations  of  components  at  u>k  a.nd  ui  where 
u>ic+0Ji  = u;m,  and  the  second  term  of  Eq.  3.5  is  a convolution  in  discrete  form.  The  coupling 
coefficient  AL{(jJk,i^t)  can  be  obtained  by  multiplying  each  side  of  Eq.  3.5  by  Yl{u}k)  and 
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Figure  3.2:  Domains  for  computing  5(wfc,a;/)aiid  B{u}i,u)k  -uji) 
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Yl(u>i)  and  by  taking  the  expected  value.  That  is, 


B*{uk,u;i) 


E[\YLiuk)YL{u,)\^] 


(3.6) 


where  B*{uk,ui)  = complex  conjugate  of  the  bispectrum  B{uk,u}i).  It  should  be  noted 
that  the  domain  of  the  bispectrum  in  this  case  is  limited  in  B,  as  indicated  in  Figure  3.2. 
The  spectral  density  function  at  the  frequency  Um  can  be  obtained  from  Eqs.  3.5  and 


3.6  as 


S{Um)  = e[\YlM\'^]+  ■ E [\YLiUk)YL{ut)\^] 


= E[\YL{u;m)\ 


|5(wjt,W()|' 


u,*+£w„.  ^ 


(3.7) 


= E[\YlM\^]+  b'{uk,Ulf  ■ SiUm) 


where 


b\ujk,i^t)  = 


|5(a;fe,wj)|' 


E[\YlMYl{ojiW]  E[\Y{umW] 


(3.8) 


The  second  term  of  Eq.  3.7  represents  the  accumulation  of  energy  densities  associated 
with  interactions  that  occur  at  the  frequency  Uk  + cj/.  Note  that  E]\YL{uk)YL{^i)\^]  in 
Eq.  3.7  is  unknown  in  advance;  hence  we  may  evaluate  E[\Y {uk)Y{u;i)\^]  and  use  it  as  an 
initial  value  in  finding  b\uk^^i)  by  iteration. 

For  the  nonlinear  components  that  occur  at  the  frequency  in  Eq.  3.6 

and  B{ukj^i)  in  Eq.  3.7  axe  replaced  by  B*{ui^Uk  - u;/)  and  B{(jJi^Uk  - cj/),  respectively. 
Thus,  we  can  write 


S{Um)  = E [lyL(^m)|"]  + E b\u;,,iJk  - OJlf  • ‘?(c^m) 


(3.9) 


where 


6' (wi,  Wife  - a;/)  = 


\B{ui,Uk  - w/)| 


E[\YL{Uk)YL{ui)\'^]  E[\Y{Um)\^] 


(3.10) 


Note  that  the  domain  of  the  bispectrum  in  this  case  is  limited  to  (B  — B,)  illustrated  in 
Figure  3.2. 
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3.2  Example  of  Application 

As  an  example  of  application  of  the  method  for  separating  the  spectrum  into  linear  and 
nonlinear  components,  we  may  analyze  the  wave  records  shown  in  Figure  3.3,  which  are 
portions  of  the  wave  data  obtained  at  various  water  depths  in  the  ARSLOE  Project  during 
storm.  A brief  description  about  the  ARSLOE  project  is  presented  in  Appendix  B. 

The  three  wave  records  shown  in  Figure  3.3  were  obtained  from  locations  where  mean 
water  depths  were  21.4,  10.1,  and  2.1  meters,  respectively.  At  Gage  No.  F710,  the  deepest 
of  the  three,  the  surface  profile  was  almost  identical  to  that  observed  in  deep  water  waves. 
As  the  water  depth  decreased  from  the  top  record  to  the  bottom,  the  skewness  of  the  surface 
elevation  increased.  At  Gage  No.  D615  with  mean  water  depth  2.1  meters,  the  wave  profile 
shows  a definite  excess  of  high  crests  and  shallow  troughs  because  of  the  limited  water 
depth,  which  is  a typical  feature  of  a non-Gaussian  random  process.  Table  3.1  summarizes 
the  time,  distance  from  the  shoreline,  water  depth,  and  statistics  computed  from  the  time 

series  of  surface  elevation  at  three  locations. 

In  order  to  identify  the  nonlinear  components  of  waves,  bispectrum  and  bicoherence 
spectrum  are  obtained  from  the  time  histories  at  Gage  Nos.  D615,  D625,  and  F710.  The 
computation  of  the  bispectrum  and  bicoherence  spectrum  using  Fast  Fourier  Transform  is 
straightforward.  However,  care  must  be  taken  in  estimating  the  statistical  error  within 
a certain  significant  level.  Statistical  stability  of  bispectrai  estimates  is  obtained  by  the 
ensemble  averaging  or  frequency  merging  technique.  The  detailed  computational  procedure 

of  the  bispectrai  estimation  is  described  in  Appendix  C. 

To  understand  the  nonlinearity  of  random  waves,  the  bicoherence  squared  spectrum  is 
used  to  identify  significant  triad  wave-wave  interactions.  Figure  3.4  shows  the  bicoherence 
spectra  at  Gage  Nos.  D615,  D625,  and  F710.  To  facilitate  visualizing  the  bicoherence  plots, 
only  the  bicoherence  squared  values  above  the  80  percent  significance  level  for  zero  true 
bicoherence  are  shown  in  the  bifrequency  octant  (see  Appendix  C for  the  significance  level 
for  zero  true  bicoherence).  The  slab  thickness  of  0.09  in  Figure  3.4  indicates  the  80  percent 
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Figure  3.3:  Portions  of  wave  records  measured  during  the  ARSLOE  project  at  locations 
F710,  D625,  and  D615  where  mean  water  depths  are  21.4, 10.1,  and  2.1  meters,  respectively. 
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Table  3.1:  Summary  of  statistics  computed  from  the  time  histories  of  surface  elevation 
obtained  at  Gage  Nos.  D615,  D625,  and  F710 


Gage  No. 

D615 

D625 

F710 

Time 

10/25/80 

10/25/80 

10/25/80 

11:35:00-11:55:00 

11:35:00-11:55:00 

11:15:00-11:35:00 

a.m. 

a.m. 

a.m. 

Distance  from 

the  shoreline 

60 

456 

12,000 

in  meters 

Mean  water  depth 

in 

2.1 

10.07 

21.43 

meters 

Moments  of  surface 

Elevation 

m2 

0.2167 

0.8507 

1.0795 

m3 

0.1042 

0.5997 

-0.0148 

7714 

0.1958 

2.5345 

3.6721 

Variance 

0.2167 

0.8507 

1.0795 

A3  (skewness) 

1.0327 

0.7643 

-0.0132 

A4(kurtosis  — 3) 

1.1713 

0.5020 

0.1511 

Max.  and  Min.  surface 

elevation  in  meters 

Vmax 

2.25 

4.04 

3.48 

ymin 

-1.02 

-2.19 

-3.61 
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significant  level  for  zero  bicoherence  squared.  The  frequency  resolution  Aw  is  27tx 0.0078  rps, 
therefore,  values  in  the  k and  / axes  of  bifrequency  plane  in  Figure  3.4  multiplied  by  Aw 
are  radian  frequencies  wjt  and  w/,  respectively.  The  value  of  minimum  density  w,,  defined 
in  Figure  3.1,  is  assumed  to  be  7Aw  rps  at  Gage  Nos.  D615,  D625,  and  F710.  At  Gage 
No.  F710  in  Figure  3.4,  two  groups  of  peaks  are  observed  in  domains  B,  and  {B  - 5,) 
as  illustrated  in  Figure  3.2,  which  indicate  sum  and  difference  wave-wave  interactions,  re- 
spectively. For  example,  a peak  located  in  the  bifrequency  plane  B,  at  frequences  w*  and 
w/  represents  a triad  sum  interaction,  wjt  -|-  w/  w^;  meaning  that  the  sum  interaction 
between  the  spectral  wave  components  at  u>k  and  ui  generates  a new  spectral  component 
at  Wm-  Furthermore,  the  bicoherence  squared  value  6^(wjt,Wi)  represents  a fraction  of  non- 
linear energy  at  Um  due  to  spectral  components  at  w^  and  w/.  Likewise,  a peak  located  in 
the  bifrequency  plane  (B  - B,)  at  frequencies  w/  and  (wfc  - W()  represents  a triad  differ- 
ence interaction,  Wfc  - w,  w„.  As  can  be  seen  in  Figure  3.4,  the  number  of  wave- wave 
interactions  increases  dramatically  as  the  water  depth  decreases. 

The  bicoherence  spectrum  can  be  used  for  discrimination  between  nonlinearly  coupled 
components  and  linear  components  and  for  evaluating  the  fraction  of  wave  energy  in  a 
measured  spectrum.  However,  it  does  not  provide  any  direct  information  on  the  relative 
contribution  of  the  interactions  to  the  total  skewness.  This  information  is  obtained  from 
the  real  part  of  the  bispectrum.  Figure  3.5  shows  the  real  part  of  the  bispectrum  for 
the  same  Gage  Nos.  as  Figure  3.4.  It  is  weU  known  that  the  volume  of  the  real  part  of 
the  bispectrum  represents  the  mean  cube  of  the  record.  The  real  part  of  the  bispectrum 
B(wfc,W|)  represents  the  contribution  to  the  mean  cube  of  the  record  from  the  product 
of  three  Fourier  components  whose  resultant  frequency  is  zero  (w*;  ->r  = 0).  For 

example.  Figure  3.5  shows  a large  peak  located  around  B(wjfe,w,)  = B(12Aw,  12Aw),  this 
peak  indicates  that  the  self-interaction  of  primary  components  at  12Aw  generates  the  first 
harmonic  frequency  at  24Aw  in  the  auto  spectrum.  In  order  to  identify  the  location  of  large 
peaks  in  the  real  part  of  the  bispectrum,  a contour  plot  of  the  bispectrum  at  Gage  No. 
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D615  is  shown  in  Figure  3.6.  The  figure  indicates  that  the  second  largest  peak  is  located  at 
= (24Au;,  12Aw).  This  interaction  indicates  that  the  wave-wave  interaction  between 
the  primary  and  the  first  harmonic  components  generates  a second  harmonic  component  at 
36Au;  in  the  auto  spectrum. 

Figures  3.7,  3.8,  and  3.9  show  the  result  of  the  separation  of  linear  and  nonlinear 
energy  spectra  obtained  at  Gage  Nos.  F710,  D625,  and  D615,  respectively.  These  figures 
show  the  ratio  of  nonlinear  energy  to  the  total  energy  increases  as  water  depth  decreases.  As 
is  shown  in  Figure  3.7,  almost  no  energy  exists  in  low  frequency  components  (frequencies 
less  than  w,  = 7Au>  = 0.343  rps).  However,  Figures  3.8  and  3.9  show  that  nonlinear 
components  of  the  spectral  density  are  pronounced  at  low  and  high  frequencies  and  that 
no  nonlinear  component  exists  in  the  neighborhood  of  the  frequency  where  the  spectrum 

peaks. 

The  results  of  this  part  of  the  study  will  be  used  for  the  development  of  a probability 
density  function  applicable  to  the  surface  elevation  of  random  waves  in  finite  water  depth. 
(See  section  4.1  Eq.  4.3.) 
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Figure  3.4:  Bicoherence  squared  spectra  at  Gage  Nos.  F710,  D625,  and  D615  with  water 
depths  of  21.4,  10.1,  2.1  meters,  respectively. 
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Figure  3.5:  Real  part  of  bispectra  at  Gage  Nos.  F710,  D625,  and  D615  with  water  depths 
of  21.4,  10.1,  2.1  meters,  respectively. 
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Figure  3.6:  Contour  plot  of  the  real  part  of  bispectrum  at  Gage  No.  D615 
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GAGE  NO.  F710  (10/25/80  11  15-35) 


Figure  3.7:  Separation  of  linear  and  nonlinear  components  of  spectrum  at  Gage  No.  F710 
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GAGE  NO.  D625  (10/25/80  11  35-55) 


Figure  3.8:  Separation  of  linear  and  nonlinear  components  of  spectrum  at  Gage  No.  D625 
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GAGE  NO.  D615  (10/25/80  11  35-55) 


Figure  3.9:  Separation  of  linear  and  nonlinear  components  of  spectrum  at  Gage  No.  D615 


CHAPTER  4 

PROBABILITY  DISTRIBUTION  OF  WAVE  PROFILES 

This  chapter  consists  of  three  sections.  In  the  first  section,  the  application  of  the  concept 
of  Kac-Siegert’s  method  to  the  probability  distribution  of  waves  in  finite  water  depth  is 
discussed.  The  second-order  interaction  coefficients  are  evaluated  in  terms  of  bispectral  and 
linear  spectral  density  functions.  Numerically  computed  probability  density  function  and 
the  histogram  constructed  from  the  wave  data  is  compared. 

In  the  second  section,  the  derivation  of  an  asymptotic  probability  density  function  in 
closed  form  is  presented  by  approximating  the  Kac-Siegert  solution  as  a function  of  a single 
random  variable  whose  parameters  are  determined  from  information  on  cumulants  k2  and 
fca  of  the  wave  record. 

The  third  section  comprises  comparisons  between  the  newly  developed  function  and 
field  data  obtained  in  the  ARSLOE  project  during  storm  for  water  depths  of  2.1,  10.07,  and 
21.43  meters. 

4.1  Application  of  the  Kac-Siegert  Method  to  Probability  Distribution  of  Nonlinear  Waves 

As  stated  in  Chapter  1,  the  concept  of  the  Kac-Siegert  method  is  applied  for  representing 
the  statistical  properties  of  nonlinear  waves.  For  this,  we  first  write  the  surface  profile  of 
second-order  random  waves  in  finite  water  depth  as  follows: 


3/(0 


N 

Re^CkC 
k=i 
N N 


i{LJk 


(4.1) 


+ Re  ^ ^ CfcC/ 
/=1 


9kl 


where  u = radian  frequency,  e = phase  lag. 
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Qki  is  a coefficient  associated  with  + w;),  which  is  the  sum  of  the  interaction  between 
two  frequency  components  u)k  and  u>i,  while  r^i  is  associated  with  the  difference  of  the 
interaction  between  Uk  and  w;.  Since  there  is  a phase  shift  between  the  two  frequencies  Uk 
and  w;,  in  general,  these  interaction  coefficients  are  complex  numbers. 

As  reviewed  in  Chapter  2,  Eq.  4.1  has  often  been  considered  for  analysis  of  nonlinear 
random  waves  in  which  the  interaction  coefficients  are  evaluated  from  the  weakly  nonlinear 
solution.  However,  the  weakly  nonlinear  solution  cannot  be  applied  to  the  present  problem. 
Instead,  we  will  obtain  the  interaction  coefficients  through  spectral  analysis.  For  this,  let 
us  write  the  first  term  of  Eq.  4.1  as  follows: 

N N 

^ Ck{cos{ukt  + (k)  + i sin{oJkt  + ffc)}  (4-2) 

k=i  k=\ 

We  assume  that  the  linear  component  t/i(f)  is  a Gaussian  random  process.  By  ignoring 
the  factor  pg,  \c\  represents  the  spectral  density  denoted  by  2SL{uk)Au,  where  the 
subscript  L stands  for  the  spectral  density  associated  with  the  linear  wave  components. 
That  is,  the  spectral  density  z,t  uj  = Uk  is  written  by 

r I ^ 

Sk  = {25L(wjt)Aw}*/^  ~ \ J j 

Then  by  defining 

Cfc  COs{Ukt  + €k)  = 

Ck  s\n{ukt  + €k)  = 

Eq.  4.2  can  be  written  as 

N 

yi{t)  = ReY^Skiuk  + ivk)  (4.5) 

fc=i 

Note  that  the  left  side  of  Eq.  4.4  is  associated  with  a Gaussian  distribution  with  zero 

mean  and  variance  s|;  hence,  Uk  and  Vk  are  statistically  independent  random  variables. 

Similarly,  the  second  term  of  Eq.  4.1  can  be  written  as 

N N 

2/2(0  = i2e  ^ ^ Sfc  5/  { qki  {uk  + i Vk)(ui  + ivi)  + rki{uk  + i Vk)(ui  - i vi)  } (4.6) 

k=i  1=1 


Sk  Uk  (4.4) 

Sk  Vk 
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where 

St  = {2Sl{uji)Auj}^  . (4.7) 

By  taking  the  real-part  of  Eq.  4.5  and  4.6,  we  have 

y[t)  = yi(0  + J/2(0  (4-8) 

N N N 

= SkUk  + ^ Y^{skSi{qki  + rki)ukUi  -f  SkSi{qki  - rkt)vkVt} 

k=l  k=l 1=1 

To  evaluate  Sk  and  5/  in  Eq.  4.8,  it  is  necessary  to  separate  the  auto  spectrum  of  waves 
into  two  parts;  one  associated  with  linear  components,  the  other  associated  with  nonlinear 
components  of  the  waves.  The  method  of  separation  of  linear  components  from  the  auto 
spectrum  is  presented  in  the  previous  chapter.  The  interaction  coefficients  q^i  and  r^i  in 
Eq.  4.8  will  be  obtained  with  the  aid  of  the  bispectrum  of  the  waves  y{t). 

Let  us  consider  the  bispectrum  The  bispectrum  is  most  commonly  evaluated 

in  the  fundamental  region,  which  is  the  octant,  the  domain  defined  by  0 < cj/  < and 
0 < u^k  ^ ^ ^ shown  in  Figure  3.2.  The  volume  under  the  bispectrum  is  equal  to  the  3rd 
moment  of  y{t)  when  the  mean  value  of  y{t)  is  zero.  That  is,  the  3rd  moment  is  given  by 

N N 

E [{l/(0}^]  Re{B{uk,iJ,)}  (4.9) 

k=i 1=1 

On  the  other  hand,  the  third  moment  of  y(t)  = y\{t)  -f-  2/2(0  can  be  obtained  as 


= 3E  [yi(0^J/2(0  + E {j/2(0}^ 


(4.10) 


Note  that  since  yi(t)  is  a normal  distribution  with  zero  mean,  that  E[{yi{t)}^]  and  E[yi{t){y2{t)}^] 
are  zero. 

The  linear  component  yi{t)  is  usually  much  greater  than  the  nonlinear  component  ^2(05 
therefore,  the  second  term  of  Eq.  4.10  may  be  neglected  in  comparison  with  the  first  term. 
Hence,  we  have  from  Eq.  4.8, 


N N 


E [{y(0}^]  = '^E  [2/i(0^1/2(0]  = + rki) 


k=l l=l 


(4.11) 
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Then,  Eqs.  4.9  and  4.11  yield 


AT  N 


N N 


X)  + ni)  = iz  n 


Re{B(oj,,,oJi)} 


k=i 1=1 


k=l 1=1 


,2,2 


(4.12) 


As  stated  in  the  previous  section,  the  domain  of  the  spectrum  applicable  for  the  inter- 
actions associated  with  the  sum  of  two  frequency  components  is  B,  shown  in  Figure  3.2. 
Therefore,  from  Eq.  4.12,  the  interaction  coefficient  qki  can  be  obtained  as 


<1M  = (4.13) 

^k^l 

Similarly,  the  interaction  coefficient  r^/  can  be  evaluated  as 

rjb/  = -^{B{ui,(jJk  - <^i)  - ■B,(wjt,u;()}  (4.14) 

4^1 

The  functional  relationship  between  the  interaction  coefficients  qki,  Tki  and  the  param- 
eters fij  and  Xj  of  Kac-Siegert’s  solution  was  derived  in  matrix  form  by  Langley  (1987). 
Although  Langley  obtained  the  interaction  coefficients  from  the  potential  theory,  the  func- 
tional relationship  between  (?*/,  Tki)  and  (/3j,  Xj)  can  stiU  be  applied  to  the  present  problem. 
Langley  presented  Eq.  4.8  in  the  following  matrix  form: 


y(t)  = ^ u + u^(Q  + R)  M.  + - R)  H 


(4.15) 


where  Q and  R are  real  symmetric  matrices  whose  kl-th.  components  are  SkSiqki  and  SkSirki, 
and  s,  u,  and  v are  vectors  whose  kl-ih  components  are  Sjt,  ujt,  and  Vk,  respectively.  We 
may  write  two  matrices  in  Eq.  4.15  as 


Q^  + R = MifAimi  (4.16) 

Q_  — R = A2W2 

where  the  column  vectors  of  Wj  {j  =1,2)  are  normalized  eigenvectors.  They  are  orthonor- 
mal  vectors  satisfying  the  condition  inj  • wj  — /,  where  / is  the  identity  matrix.  The 
elements  of  diagonal  matrices  Aj  (y  = 1,2)  are  eigenvalues  of  matrix  Q_-\-  R and  — R, 


33 


respectively.  Since  ^ + R and  Q - are  iV  by  iV  symmetric  matrices,  there  are  a total  of 
2N  eigenvalues. 

Futhermore,  by  writing 

f j = 1,2,- --iV 

Z,-  = { (4.17) 

I {wU)j-N.  j = {N  + l)r-2N 

where  Zj  is  an  independent  Gaussian  random  variable  with  zero  mean  and  unit  variance, 
Langley  has  shown  that  Eq.  4.15  can  be  reduced  to  the  Kac-Siegert’s  formulation  given  in 
Appendix  A Eq.  A. 2 with  the  following  relationships: 


(4.18) 

1 0’ 

j = {N  + l),---2N 

f (Ai)i 

j = l,2,---N 

A,= 

(4.19) 

[ {^2)31 

j = (iV  + l),---2iV 

By  employing  the  parameters  (3j  and  Xj  obtained  in  Eqs.  4.18  and  4.19,  the  probability 
density  function  of  y(t)  can  be  derived  numerically  from  Eq.  A. 2.  Hereafter,  y{t)  will  be 
written  as  a random  variable  Y,  The  most  convenient  way  to  derive  the  probability  density 
function  is  through  the  characteristic  function  as  suggested  by  Kac  and  Siegert.  That  is, 
the  characteristic  function  of  Eq.  A. 2 is  given  by 


2N 


2N 


2*2 


m = jQ(i  - g 

and  therefrom  the  probability  density  function  of  Y can  be  numerically  derived  by 


(4.20) 


!(y)  = ^ r *(<)  * 


(4.21) 


The  method  stated  above  is  applied  to  the  data  obtained  at  Gage  No.  D615  shown  in 
Figure  3.3  in  which  fJj  and  Xj  are  evaluated  by  Eqs.  4.18  and  4.19  for  N = 128  points.  A 
comparison  between  the  probability  density  function  and  the  histogram  constructed  from 
the  data  is  shown  in  Figure  4.1.  Excellent  agreement  between  them  can  be  seen  in  the 

figure. 


PROBABILITY  DENSITY  IN  1/METER 
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Figure  4.1:  Compaiison  of  probability  density  function  obtained  fjom  Kac-Siegert 
and  histogram  constructed  from  data 
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4.2  Derivation  of  Asymptotic  Probability  Density  Function 

As  demonstrated  in  the  previous  section,  Kac-Siegert’s  solution  satisfactorily  yields  a 
probability  density  function  applicable  to  nonlinear  random  waves  in  finite  water  depth.  The 
probability  density  function,  however,  is  not  given  in  closed  form.  One  way  to  derive  the 
probability  density  function  in  closed  form  is  to  present  Kac-Siegert’s  solution  (Eq.  A. 2)  as 
a function  of  a single  random  variable  instead  of  the  summation  of  the  standardized  normal 
distribution  and  its  squared  quantity.  For  this,  let  us  present  Eq.  A. 2 as 

Y = U + aU^  (4.22) 

where  is  a constant  (unknown)  and  (7  is  a normal  variate  with  mean  and  variance 
al,  both  of  which  are  also  unknown.  The  value  of  these  unknowns  will  be  determined  by 
equating  the  characteristic  function  of  Eq.  4.22  with  that  of  Kac-Siegert’s  solution  given  in 

Eq.  4.20. 

In  order  to  obtain  the  characteristic  function  of  Eq.  4.22,  let  us  write 


(4.23) 


where 

Uo  = y/a  ^ . (4.24) 

Since  is  a normal  variate  with  mean  and  variance  cr^,  Uo  follows  the  normal  distribution 
with  mean  ^/aa^  + (l/(2y^)  sind  variance  aa^.  Therefore,  the  characteristic  function  of 

Uq  becomes 


4>{t)  = 


1 


y/1  — i2aalt 


exp 


it  (yan,  + 

1 — i2aa‘it 


The  characteristic  function  of  the  second  term  of  Eq.  4.23  is  given  by 


<f>{t)  — exp 


(4.25) 


(4.26) 
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Thus,  the  characteristic  function  of  Y is  the  product  of  Eqs.  4.25  and  4.26  which  after  some 
mathematical  manipulation  can  be  written  as  follows: 


<f>{t)  = 


1 


— i2aalt 


exp  < ( a/i*  + /X.  + 


it 


1 — i2aalt 


(4.27) 


By  taking  the  logarithm  of  the  characteristic  function  and  by  expanding  it  in  a series 


of  (it),  we  have 


ln(f>{t)  = 


+ 


{aal  + aul  + /x,)(it)  + ja^cr^  + y2a(r^(a^^  + ^.)  \ (it)' 

+ 4a^<T^(a/x^  + /x.)|  (it)^  + • • • 


(4.28) 


On  the  other  hand,  by  taking  the  logarithm  of  the  characteristic  function  of  Kac- 
Siegert’s  solution  (Eq.  A.2),  and  by  expanding  it  in  a series  of  (ti),  we  have 


(2N  \ /l  2N  \ 


2N 


2N 


(4.29) 


,i=i 


i=i 


Note  that  the  first  term  is  = 0»  since  the  mean  value  of  Y is  zero.  Comparing  Eqs.  4.28 
and  4.29,  the  following  relationship  can  be  derived: 


o<t2  + aul  + ^,  = 0 


2N  2N 

^ i=i  i=i 


— ^2 


,2_4 


at  - 2aV; 


(4.30) 


2N  2N 

3 E + 4 E ^-  = - Sa^at 

j=i  i=i 

By  using  the  values  of  /3j  and  Xj  evaluated  in  the  previous  section,  we  can  determine 
the  unknown  parameters  a,  and  from  Eq.  4.30.  Since  U in  Eq.  4.22  is  now  a normal 
variate  with  known  mean  and  variance,  the  probability  density  function  of  Y can  be  derived 


by  applying  the  technique  of  change  of  random  variables  from  U to  Y . The  density  function 


can  be  obtained  as 
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m = 

V2iracr,  ^ 


A compajison  of  the  asymptotic  probability  density  function  shown  in  Eq.  4.31  and  that 
obtained  numerically  from  the  characteristic  function  of  Kac-Siegert’s  solution  is  shown  in 
Figure  4.2  along  with  the  histogram  of  a random  process.  As  can  be  seen  in  the  figure,  the 
asymptotic  probability  density  function  agrees  well  with  that  obtained  from  the  Kac-Siegert 
solution.  Unfortunately,  however,  the  density  function  vanishes  at  the  point  y = -(l/4a) 
due  to  a singular  point  involved  in  the  density  function.  Because  of  this  drawback,  Eq.  4.31 
cannot  be  considered  as  an  appropriate  density  function. 

Next,  let  us  present  the  functional  relationship  between  Y and  U given  in  Eq.  4.22 
inversely  such  that  the  random  variable  U is  expressed  as  a function  of  F as  follows: 

U=  — (l-  (4.32) 

7a  v ' 

where  7 is  a constant;  1.28  for  y > 0 and  3.00  for  y < 0.  Justification  for  selecting  these 
constant  values  of  7 is  given  in  Appendix  D.  It  may  suffice  here  to  say  that  these  constants 
are  determined  from  dimensionless  quantities  of  Y and  U,  and  that  the  constants  are  valid 
even  for  waves  with  very  strong  nonlinear  characteristics. 

It  is  noted  that  the  values  of  7 are  different  for  positive  and  negative  y-values.  This 
results  in  a slight  difference  in  the  slope  of  the  probability  density  function  at  y = 0,  although 
the  density  function  is  continuous  at  this  point. 

Figure  4.3  shows  a comparison  of  Eq.  4.22  and  Eq.  4.32  for  a = 0.326.  The  a-value  is 
obtained  from  Eq.  4.30  for  the  non-Gaussian  random  process  shown  in  Figure  3.3,  Gage  No. 
D615  which  has  a very  large  skewness,  1.03.  The  random  process  may  be  considered  as  if  it 
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Figure  4.2:  Comparison  between  Kac-Siegert  and  asymptotic  probability  density  functions 
and  histogram  constructed  from  data 
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were  random  waves  with  very  severe  nonlinearity.  As  can  be  seen  in  the  figure,  Eq.  4.32  well 
represents  Eq.  4.22  except  in  the  negative  range  of  y where  Eq.  4.22  has  a singular  point. 

By  using  the  functional  relationship  given  in  Eq.  4.32,  we  may  derive  the  probability 
density  function  of  Y from  the  random  variable  U which  obeys  the  normal  distribution  with 
mean  /z,  and  variance  a^.  The  change  of  random  variable  technique  yields  the  probability 
density  function  of  T as  follows: 


where 


1.28 


J/  > 0 


7 = 


\ 3.1 


(4.34) 

.00  3/  < 0. 

By  applying  the  method  for  determining  the  limit  of  an  indefinite  function,  it  can  be 
easily  proved  that  Eq.  4.33  reduces  to  a normal  probability  density  function  with  mean  /z. 
and  variance  cr^  if  a = 0;  namely,  for  linear  random  waves. 

It  was  stated  earlier  that  a,  /z,,  and  cr^  are  determined  from  Eq.  4.30  using  j3j  and  \j 
values  of  Kac-Siegert’s  solution.  However,  this  method  of  evaluating  a,  /z,,  and  al  can  be 
greatly  simplified  if  we  consider  the  cumulants  of  the  distribution.  That  is,  the  cumulant 
generating  function,  denoted  by  rp{t),  of  the  Kac-Siegert  solution  can  be  obtained  from 
Eq.  4.20  as  follows: 


1 3t^  1 

xj){t)  = ln4>{t)  = -- 2 l-t'2A 

^ i=i  i=i  ^ 

By  expanding  the  first  and  second  terms  in  a series  of  (i^),  we  have 


t 


(4.35) 


(2N  \ 1/2^  2N  \ 

(“)  + 2!  ( 1 


+ 


, / 2N  2N  \ 


(4.36) 


On  the  other  hand,  the  cumulant  generating  function  can  be  written  in  terms  of  cumu- 


40 


Figure  4.3:  Comparison  of  F = U + aU'^  and  t/  = (1  - exp{-a7F})/7o  for  a - 0.326 
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lants  as  follows: 


(4.37) 


where  ks  = cumulants. 


From  comparison  of  Eqs.  4.36  and  4.37,  we  have  the  following  relationships: 


i=i 

2N  2N 


2N 


^2  = 13  ^ 13 


(4.38) 


j=i  i=i 

2N  2iV 


2N 


^3  = 6 ^ /3j  Aj  + 8 13 


Therefore,  Eq.  4.30  can  be  written  as 


aal  + + M*  = ^ 

< al  - 2a^a^  - k2 


(4.39) 


2a<r;J(3  — 8a*<rJ)  = 


Since  the  mean  value  of  the  surface  elevation  Y is  zero  for  the  present  problem,  we 
have  ki  = 0,  and  thereby  fcj  and  k^  are  equal  to  the  sample  moments  E[y^\  and  E[y% 
respectively.  Thus,  we  can  determine  the  unknown  parameters  a,  /i*,  and  cr,  by  simply 

evaluating  the  sample  moments  from  the  wave  profile  record  y{t). 

4.3  AppUcation  of  the  ProbabiUty  Density  Function  to  Field  Data 

The  newly  developed  probability  density  function  applicable  to  surface  profiles  of  ran- 
dom waves  in  finite  water  depth  is  compared  with  the  histograms  constructed  from  field 
data  obtained  at  various  water  depths.  The  wave  data  used  are  shown  in  Figure  3.3. 

As  mentioned  in  the  previous  section,  the  parameters  a,  /x*,  and  al  are  determined 
from  Eq.  4.39  with  the  cumulants  k2  and  ks  computed  from  the  record.  Table  4.1  shows 

the  parameters  computed  at  three  different  water  depths. 

Figures  4.4,  4.5,  and  4.6  show  comparisons  of  the  probability  density  function  given 

in  Eq.  4.33  and  the  histogram  constructed  from  wave  record  at  Gage  Nos.  D615,  D625,  and 
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Table  4.1:  Parameters  computed  from  the  cumulants  of  surface  elevation  at  Gage  Nos. 
D615,  D625,  and  F710 


Gage  No. 

D615 

\ 

D625 

F710 

Cumulants  of  surface 
Elevation  in  meters 

k2 

0.2167 

0.8507 

1.0795 

ks 

0.1042 

0.5997 

-0.0148 

Parameters 

Computed 

a 

0.3716 

0.1351 

-0.0021 

-0.089 

-0.121 

0.002 

0.2315 

0.8790 

1.0795 

F710,  respectively.  Included  also  in  the  figures  are  the  normal  probability  density  function 
with  zero  mean  and  variance  evaluated  from  the  records.  Figure  4.4  shows  a comparison  for 
waves  in  shallow  water  areas  (mean  water  depth  2.1  meters).  As  can  be  seen,  the  histogram 
deviates  from  the  normal  density  function  to  a great  extent;  but  the  agreement  between  the 
histogram  and  the  newly  developed  probability  density  function  is  excellent. 

Figure  4.5  shows  the  case  for  waves  in  water  depth  10.1  m in  which  the  skewness 
representing  the  degree  of  nonlinearity  of  waves  is  less  than  that  of  Figure  4.4.  Again,  the 
probability  density  function  shows  excellent  agreement  with  the  histogram. 

Figure  4.6  shows  a comparison  for  a water  depth  of  21.4  m.  In  this  case,  the  distribution 
of  surface  elevation  is  the  normal  distribution  as  shown  in  the  figure  because  the  water  depth 
is  relatively  deep.  The  histogram  and  the  newly  developed  probability  density  function  agree 
very  well  in  this  case. 


PROBABIUTY  DENSITY  IN  l/M 
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GAGE  NO.  D615  (10/15/80  11  35-55) 


Figure  4.4;  Comparison  of  newly  developed  probability  density  function  and  histogram 
constructed  from  data  at  Gage  No.  D615 


PROBABIUTY  DENSITY  IN  l/M 
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GAGE  NO.  D625  (10/25/80  11  35-55) 


Figure  4.5:  Comparison  of  newly  developed  probability  density  function  and  histogram 
constructed  from  data  at  Gage  NO.  D625 


PROBABIUTY  DENSITY  IN  1/M 
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GAGE  NO.  F710  (10/25/80  11  15-35) 


Figure  4.6:  Comparison  of  newly  developed  probability  density  function  and  histogram 
constructed  from  data  at  Gage  No.  F710 


CHAPTER  5 
CONCLUSIONS 


This  paper  presents  the  results  of  a study  on  nonlinear  analysis  of  random  waves  in  finite 
water  depth,  in  particular,  this  study  seeks  to  (1)  clarify  nonlinearities  in  wave  spectra,  and 
(2)  develop  a probability  density  function  in  closed  form  which  is  applicable  to  random 
waves  in  finite  water  depth. 

In  the  first  phase  of  the  study,  a method  to  separate  the  second-order  spectrum  from  the 
measured  spectrum  was  developed.  The  method  employed  bicoherence  spectrum  through 
which  nonlinear  characteristics  of  waves  in  finite  water  depth  were  clarified.  The  method 
was  applied  to  field  data  at  different  water  depths  obtained  in  the  ARSLOE  Project  during 

storm. 

From  the  results  of  a separation  of  the  second-order  spectrum  from  the  measured  spec- 
trum, it  was  found  that  no  nonlinear  components  existed  near  the  spectrum  peak  frequency, 
Wp.  In  shallow  water  depth,  nonlinear  components  of  the  spectral  density  were  pronounced 
at  high  and  low  frequencies  which  indicated  that  a considerable  number  of  occurrences  of 
wave- wave  interactions  were  associated  with  various  pairs  of  two  frequency  components. 

It  was  also  found  that  the  ratio  of  nonlinear  energy  to  the  total  energy  increased  with 
decreasing  water  depth.  This  nonlinear  energy  was  associated  with  mainly  self-interactions 
at  frequencies  near  ojp.  The  interactions  between  frequencies  in  the  neighborhood  of  Wp  and 
their  harmonic  components  (2u;p  or  3u>p)  became  pronounced  with  decreasing  water  depth. 

In  the  second  phase  of  the  study,  the  probability  density  function  applicable  to  random 
waves  in  finite  water  depth  was  developed  in  closed  form  by  applying  the  concept  of  Kac- 
Siegert’s  method  for  evaluating  the  response  of  a nonlinear  mechanical  system.  That  is, 
the  Kac-Siegert  formula  was  expressed  in  terms  of  a random  variable  which  obeys  a normal 
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distribution  with  parameters  evaluated  from  information  on  cumulants  of  the  wave  record. 
Then,  by  applying  the  transformation  technique  of  random  variables,  the  probability  density 
function  applicable  to  nonlinear  random  waves  was  developed  in  closed  form.  Comparisons 
between  the  newly  developed  probability  density  function  and  the  histograms  constructed 
from  wave  records  at  various  water  depths  obtained  in  the  ARSLOE  Project  during  strom 
showed  excellent  agreement. 


APPENDIX  A 
KAC-SIEGERT’S  METHOD 


The  response  of  a nonlinear  system  is  often  presented  by  Volterra’s  stochastic  series 
expansion  given  by 


yi{t)  + y2{t) 


(A.1) 


j h\{T)x{t  — r)dT  + j J h2{ri,T2)x{t  — Ti)x{t  — T2)d,Ti  dr2  + • • • 


where 


yi{t),y2{t) 

hi{r) 

h-2{Ti,T2) 


linear  and  nonlinear  component  of  the  response,  respectively, 
impulse  response  function, 
second- order  impulse  response  function 


For  this  expansion,  Kac  and  Siegert  (1947)  have  shown  that  the  response  y{t)  can  be 
obtained  as  the  sum  of  the  standardized  normal  variates  and  its  squares  as  follows; 


y(t)  = Z,-  + A,Z?)  (A.2) 

j 

Here,  the  parameters  /3j  and  Xj  are  evaluated  by  finding  the  eigenfunction  and  eigenvalues 
of  the  integral  equation  given  by 

J K{ui,U2)  du}2  = Aj  (■^•3) 

where 

K{uJi,UJ2)  = H{uJi,UJ2)\/S{ui)  S{U2) 

H{u}i,u>2)  = second-order  frequency  response  function 

5(a;)  = input  spectral  density  function 
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V>j(w) 


orthogonal  eigenfunction  = 


/ =1  for  j = A: 

0 otherwise 


Since  y{t)  is  represented  in  terms  of  a standardized  normal  random  variable,  the  proba- 
bility density  function  of  y{t)  can  be  derived  through  the  characteristic  function  of  Eq.  A. 2. 
This  method  developed  by  Kac  and  Siegert  is  a novel  and  meticulous  approach  for  deriv- 
ing the  probability  density  function  of  a second-order  nonlinear  response.  Unfortunately, 
however,  the  probability  density  function  cannot  be  obtained  in  closed  form.  It  is  also  un- 
derstood that  the  parameters  fij  and  Xj  in  Eq.  A. 2 are  evaluated  from  information  of  either 


the  second-order  impulse  response  function,  h(tiyt2)  or  the  second-order  frequency  response 
function,  H{ui^U2). 


APPENDIX  B 
ARSLOE  WAVE  DATA 


The  wave  data  used  in  the  present  study  were  obtained  by  the  Coastal  Engineering 
Research  Center  (CERC),  U.  S.  Army,  at  its  Field  Research  Facility,  located  at  Duck,  North 
Carolina.  The  Field  Research  Facility  has  a 550-meter-long  research  pier  extending  into  the 
Atlantic  Ocean,  equipped  with  seven  Baylor  resistance- type  wave  gages  (see  Figure  B.l). 
Extensive  wave  measurements  were  made  at  this  site  for  two  months  in  1980  under  the 
Atlantic  Ocean  Remote  Sensing  Land-Ocean  Experiment  (ARSLOE)  project. 

During  October  23  - 25,  an  extratropical  cyclone  moved  directly  through  the  experi- 
mental area  with  wind  speeds  on  the  order  of  10  — 15  meters/sec.  Continuous  wave  records 
were  taken  by  the  CERC  Field  Research  Facility  during  the  storm.  Statistical  analysis  of 
these  records  were  carried  out  by  Ochi  et  al.  (1982)  at  the  Cocistal  and  Oceanographic  En- 
gineering Department,  University  of  Florida.  The  wave  data  analyzed  in  the  present  study 
are  those  recorded  at  two  locations  D615  and  D625  along  the  pier  and  one  buoy  (wave 
rider)  F710  located  offshore  along  the  extended  bne  of  the  pier.  The  location  of  the  wave 
gages  and  wave  riders  as  well  as  the  average  water  depth  measured  on  October  21  and  27 
by  CERC  are  given  in  Figure  B.l. 
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Figure  B.l:  Coastal  Engineering  Research  Center  (CERC)  Field  Research  Facility,  Duck, 
North  Carolina  (By  courtesy  of  CERC) 


AVERAGE  WATER  DEPTH  IN  METERS 


APPENDIX  C 

COMPUTATION  OF  THE  BISPECTRUM 


This  appendix  outlines  computational  procedures  of  the  bispectral  estimation.  Suppose 
that  the  water  surface  elevation,  y(t),  is  a sample  of  size  N. 

(1)  Divide  the  total  number  of  data  points  N into  K sets  of  segments  such  that  each  segment 
has  M data  points. 

y{^'>  = k=l,...,K  and  (C.l) 

(2)  Subtract  the  mean  value  from  each  segment. 

(3)  Apply  an  appropriate  window  function  Wj  to  each  segment  to  reduce  leakage. 

vj  = wjVj,  j = l,...,M  (C.2) 

(4)  Compute  the  complex  Fourier  coefficients  by  using  FFT  technique. 

i i = 0,...,^-l  (C.3) 

One  should  scale  Vj  to  obtain  the  correct  magnitudes  (Bendat  and  Piersol,  1980). 

=J=—  i = (C.4) 

(5)  Estimate  the  bispectrum  and  the  bicoherence  spectrum  by  averaging  bispectral  values 
from  each  segment  (ensemble  averaging). 

S(u,,c,)  = (C.5) 

i=l 

where  “ B ” denotes  an  estimator  and  the  bispectrum  estimate  obtained  by  arithmetic 
ensemble  averaging  is  an  unbiased,  minimum  variance  estimator  (Kim  and  Power,  1979). 
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For  step  1 in  the  procedure,  to  achieve  the  smallest  spectral  variance  from  a fixed  number 
of  available  sample  data  points,  we  divided  the  total  number  of  data  points  N = 2176 
(sampling  frequency  is  2 Hz)  into  segments  of  M = 256  points  with  50  percent  overlap  of 
each  segment.  Therefore,  the  number  of  segments,  K,  is 


(C.6) 


According  to  Welch  (1967),  the  50  percent  overlap  is  nearly  optimal  to  have  the  smallest 
variance.  The  overall  reduction  in  variance  is  a factor  of  11/(9 A)  compared  to  the  variance 
of  each  realization,  and  the  equivalent  degrees  of  freedom  (d.o.f.)  of  approximating  chi- 
square  distribution  is 


d.o.f.  « 


3.584 


N 

M + l 


= 3.584  X 


2176 

256-1-1 


(C.7) 


According  to  Elgar  and  Guza  (1988), the  80  percent  significance  level  for  zero  bicoher- 
ence for  30  degrees  of  freedom  is  approximaly  0.3.  Therefore,  bicoherence  squared  values 
greater  than  0.09  are  selected  for  the  separation  of  second-order  spectrum  from  the  measured 
spectrum  as  explained  in  Chapter  3.2. 


APPENDIX  D 

INVERSE  PRESENTATION  OF  Y = U + aU^ 


The  functional  relationship  between  two  random  variables  Y and  U is  given  in  Eq.  4.22 
as  follows: 

Y = U + aU^  (D.l) 

where  “a”  is  a constant  and  17  is  a normal  Vciriable  with  mean  /x,  and  variance  cr^.  The 
probability  density  function  of  Y obtained  by  using  the  functional  relationship  vanished  at 
a point  Y = -l/(4a)  due  to  a singularity  involved  in  the  density  function  of  Y.  One  way 
to  solve  this  difficulty  is  to  derive  the  density  function  of  Y by  presenting  the  functional 
relationship  given  in  Eq.  D.l  inversely  such  that  the  random  variable  U is  expressed  as  a 
function  of  Y.  For  this,  let  us  write  Eq.  D.l  in  dimensionless  form.  That  is 

W = v + V'^  (D.2) 


where  W = aY  and  V — all . 

Then,  we  may  write  V as  a function  of  IF  as  follows: 

V^-(l-  e-^^)  (D.3) 

7 

where  7 = constant. 

Prior  to  determining  the  value  of  7,  let  us  find  the  domain  of  W and  V appropriate  for 
non-Gaussian  random  processes  in  practice.  We  consider  first  the  positive  domain  of  V and 
W.  As  an  example,  the  random  process  shown  in  Figure  3.3,  Gage  No.  D615  demonstrates  a 
typical  feature  of  non-Gaussian  characteristics  — a definite  excess  of  high  crests  and  shallow 
troughs.  The  skewness  of  this  example  is  1.03  which  indicates  a considerable  deviation  in  the 
shape  of  the  proability  density  function  from  that  of  a Gaussian  process.  In  this  example. 
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the  maximum  positive  deviation  from  the  mean  value  is  less  than  2.1,  and  “a”  is  0.37.  Thus, 
the  range  of  W from  0 to  0.78  (V  from  0 to  0.51)  appears  to  be  sufficient  for  this  example 
of  a random  process  having  strong  nonlinearity. 

As  another  example,  let  us  consider  the  nonlinear  turbulent  wind-  induced  drag  force 
on  a structure.  The  instantaneous  wind  speed  at  height  z above  the  sea  surface  can  be 
written  as 

U,  = U,  + w(t)  (D.4) 


where  Uz  = mean  wind  speed  at  height  z 
w{t)  = turbulent  wind  speed. 

The  wind-induced  drag  force,  Fq,  acting  on  the  structure  with  a projected  area  of  A 
can  be  evaluated  by 


pCDA{U  + w{t)}^ 


^ w{t)  -b  w{t)  |u;(t)|  = a in  -b  5 tn  |u;| 


rr 


(D.5) 


where 


p = air  density 
Cd  = drag  coefficient 

1 — 2 

Fd  = drag  force  due  to  mean  wind  speed  U = - p Cd  AU 

a = 2Fd/U,  b=FDlU^ 

By  multiplying  by  (b/a^)  on  both  sides  of  Eq.  D.5,  we  can  present  the  drag  force  in  the 

same  dimensionless  formula  as  given  in  Eq.  D.2  with  W = FD^bfa^)  and  V = (bfa)w. 

We  may  consider  an  extremely  severe  wind  having  a speed  of  38.5  m/sec  (75  knots) 
with  a maximum  turbulent  speed  of  30.8  m/sec  (60  knots)  blowing  over  a structure  with 
a projected  area  of  2 X 10^  m^.  This  results  in  a = 107  X 10^  kg/m,  and  thereby  we  have 
a maximum  V of  0.38.  In  other  words,  the  range  of  V is  from  0 to  0.35  and  thereby  W 
ranges  from  0 to  0.47  yielding  a much  smaller  domain  of  V and  W than  that  obtained  for 
the  random  process  shown  in  Figure  3.3,  Gage  No.  D615. 
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Next,  let  us  choose  the  7—  value  in  Eq.  D.3  which  yields  the  closest  agreement  with 
Eq.  D.l  in  the  positive  domain  of  W and  V discussed  above.  Results  of  computation  show 
that  7 = 1.28  in  Eq.  D.3  as  shown  in  Figure  D.l  in  the  specified  domain. 

The  negative  domain  of  W and  V is  obtained  through  a method  similar  to  that  consid- 
ered for  the  positive  domain.  For  the  non-  Gaussian  random  processes  shown  in  Figure  3.3 
Gage  No.  D615,  the  smallest  negative  deviation  from  the  mean  value  is  less  than  —1.0, 
and  thereby  the  range  of  W is  from  0 and  -0.37.  Note  that  the  probability  density  func- 
tion of  V obtained  by  using  Eq.  D.l  has  a singular  point  at  y = -l/(4a)  and  hence  at 
]V  = —ay  = —0.25.  Therefore,  it  is  appropriate  to  choose  a 7— value  in  Eq.  D.3  which  best 
fits  Eq.  D.l  for  W less  than  —0.25.  The  results  of  computation  indicate  that  7 = 3.00  in 
Eq.  D.3  as  shown  in  Figure  D.2  best  represents  Eq.  D.l  in  the  specified  range  of  W. 
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Figure  D.l:  Comparison  between  W = V and  — (1  — exp{— 7W^})/7  with 

7 = 1.28,  (K,  W > 0) 


Figure  D.2:  Comparison  between  W — V + and  _ (1  — exp{— 7lF^})/7  with 

7 = 3.00,  (K,  IF  < 0) 
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